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[57] ABSTRACT 

An improved method based on a simple imaging operation 
with a simple one-to-one mapping between resolution ele- 
ments of a spatial phase modulator and resolution elements 
of the generated intensity pattern is provided. According to 
the invention a method is provided for synthesizing an 
intensity pattern with low loss of electromagnetic energy, 
spatial modulation of electromagnetic radiation with a spa- 
tial phase mask for modulation of the phase of the incident 
electromagnetic radiation by phasor values of individual 
resolution elements of the spatial phase mask, each phasor 
value being determined in such a way that the values of the 
Fourier transformed phasors attain predetermined values for 
predetermined spatial frequencies, and the phasor value of 
the specific resolution element of the spatial phasor mask 
corresponds to a distinct intensity level of the image of the 
resolution element in the intensity pattern, and a spatial 
phase filter for phase shifting of a part of the electromagnetic 
radiation, in combination with an imaging system for gen- 
eration of the intensity pattern by interference in the image 
plane of the imaging system between the part of the elec- 
tromagnetic radiation that has been phase shifted by the 
phase filter and the remaining part of the electromagnetic 
radiation. 

65 Claims, 6 Drawing Sheets 
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PHASE CONTRAST IMAGING 
FIELD OF THE INVENTION 

The invention relates to a method and a system for 
synthesizing a prescribed intensity pattern based on phase 
contrast imaging. 

BACKGROUND OF THE INVENTION 

It is well known to form an image on an illuminated 
surface of a body by absorption or blocking of energy of an 
illuminating beam. For example in an overhead projector, an 
over-head transparent absorbs or blocks part of the light 
beam of the projector whereby a large image of an overhead 
is formed on a screen. However, this results in a loss of light 
intensity as part of the emitted light from an image forming 
system is reflected or absorbed. 

To avoid loss of energy causing, e.g. loss of light intensity 
of the synthesized intensity pattern, power dissipation gen- 
erating heat in components of the system, etc., methods and 
systems have been developed wherein the phase of a light 
beam is modulated instead of the amplitude or intensity of 
the light beam, as modulation of the phase of the light beam 
do not lead to loss of energy. The phase modulation is 
followed by a conversion of the phase modulation into an 
amplitude or intensity modulation. 

A diffractive optical element, such as a holographic opti- 
cal element, may be used to generate a phase modulation. 
Then, the resulting intensity modulation at each point of a 
picture formed by conversion of the phase modulation into 
intensity modulation will depend upon the phase modulation 
values at each point of the diffractive optical element as the 
light intensity at each point of the picture is formed by a 
coherent superposition of light received from the entire 
surface of the diffractive optical element. Diffractive optical 
elements are rather complex to design for synthesis of a 
prescribed intensity pattern. 

Imaging methods and systems may also be used in 
connection with phase modulation. These methods and 
systems are characterized by the fact that the intensity of a 
point of a picture formed by conversion of phase modulation 
into intensity modulation will depend upon the phase modu- 
lation value of one point of the phase modulator only as this 
point is imaged onto the picture point in question by the 
imaging system. This one-to-one relationship makes the 
design of phase modulators in these systems simple. Meth- 
ods and systems of this kind are named phase contrast 
imaging methods and systems. 

Phase contrast imaging methods were originally devel- 
oped within the field of microscopy. Many objects of interest 
in microscopy are largely transparent, thus absorbing little or 
no light. When light passes through such an object, the 
predominant effect is the generation of a spatially varying 
phase shift which can not be seen by a human as the eye of 
a human responds to light intensity and colour and does not 
respond to the phase of light. 

In 1935, Fritz Zernik proposed a phase contrast technique 
which rests on spatial-filtering principles and has the advan- 
tage that the observed intensity is linearly related to the 
phase shift introduced by the object. 

Suppose that a transparent object with amplitude trans- 
mittance 

is coherently illuminated in an image-forming system. For 
simplicity, a magnification of unity is assumed and the finite 
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extent of the exit and entrance pupils of the system is 
neglected. Further, a necessary condition to achieve linearity 
between phase shift and intensity is that the phase shift § be 
less than 1 radian, in which case the amplitude transmittance 
s can be approximated by 

t&yhUMw) (2) 

The terms of order $ 2 and higher are neglected in this 
approximation. It is seen that the first term of (2) leads to a 
strong wave component that passes through the sample 
without change, while the second term generates weaker 
diffracted light that is deflected away from the axis of the 
system. 

The image produced by a conventional microscope can be 
written 

/-|1 + ><M 2 -1 (3) 

where the term (j> 2 has been approximated by zero. It is seen 

20 that the diffracted light is not observable because it is in 
phase quadrature with the strong background. As Zernik 
recognized that the background is brought to a focus on-axis 
in the focal plane while the diffracted light — containing 
higher spatial frequencies — is spread away from the focal 

25 point, he proposed that a phase -changing plate be inserted in 
the focal plane to modify the phase relation between focused 
and diffracted light. 

The phase-changing plate can consist of a glass substrate 
on which a small transparent dielectric dot has been coated. 

30 The dot is placed at the center of the focal plane and has a 
thickness and index of refraction such that it retards the 
phase of the focused light by either n/2 radians or 3x/2 
radians relative to the phase retardation of the diffracted 
light. In the former case the intensity in the image plane 

35 becomes 

/=|«pD-(rt/2)+-y+l 2 -l/Yi+<t>)l 2 -i+2* (<*) 

while in the latter case 

40 /-|«ii/(3j^)+;<i>l 2 -l/(-i+4>)l a -i-24> ( 5 > 

Thus, the image intensity has become linearly related to 
the phase shift ((>, When the phase of the background is 
retarded by jt/2, the result is known as positive phase 

45 contrast, while a 3ji/2 retardation is said to yield negative 
phase contrast. 

It is seen that the method described above leads to a phase 
contrast imaging method that provides a small phase signal 
that is superimposed on a large DC-component. This leads 

50 to an important disadvantage of the method because, 
typically, it will be necessary to attenuate the DC-component 
to enhance the information contained in the phase modulated 
signal. However, the attenuation of the DC-component leads 
to loss of energy. This kind of filtering is usually denoted 

55 Dark Field Filtering. 

It is another disadvantage of the phase contrast imaging 
method described above that it is based on the assumption 
that the phase shift $ is less than 1 radian which is very often 
not fulfilled in practical real-life applications. However, the 

60 theory is still applied to such applications, disregarding the 
fact that the basic assumption is not fulfilled, and this leads 
to non -optimized technical solutions. 

In EP 0 657 760 a phase contrast imaging system is 
disclosed in which an image simulation and projection 

65 system is based on the Texas Instrument flexure beam digital 
mirror device (DMD). The flexure beam DMD is used for 
analog phase modulation of reflected light and the phase 
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modulation is converted to amplitude modulation utilizing a 

phase contrast imaging method. The flexure beam DMD ^ y f ^ Xt y f ^ y f (6) 

provides a flicker-free modulated wave and accordingly, y ^ =e ~ 1 ^ j\ 1 + ^ + ■•• 
optical image sensor synchronization is not needed. The 

system disclosed operates according to the Zernike method s . 4 , iL . ... iU iL . 

' A . . . , , L j . . . - , into account and, thus, is not based on the assumption that 

and, thus includes the corresponding disadvantages ^ phage shift ^ fc ^ ^ x fadian , t fc important to flote 

descnbed above. that each term of the Tay i or > s ser j es contribute to the 

Another example of a phase contrast imaging system is DC-value of the function t(x,y). This fact is not recognized 

disclosed in GB 2 199 716, wherein an optical guide-beam 1Q in the art as the DC-value has until now been believed to be 

projector for a missile guidance system is disclosed that represented by numeral 1 in equation (2). 

provides a spatially intensity modulated guide-beam. A It is still another object of the present invention to provide 

spatial phase modulator is used to generate the guide-beam. an improved method and apparatus for phase contrast imag- 

The phase encoding of the spatial phase modulator consti- ing without the need of attenuating the DC-component of the 

tutes a periodic square- wave modulation (50% duty cycle) 35 signal to enhance the information contained in the phase 

of two phase values 0 and ji/2. The phase modulation is modulated signal. 

converted into an amplitude modulation by Fourier trans- 11 is yet another object of the present invention to provide 

forming lenses and a phase plate providing a phase shift of an improved method based on a simple imaging operation 

the background signal by n/2- A method for synthesizing the with a simple one-to-one mapping between resolution ele- 

specific intensity pattern of the optical guide-beam based on 20 ments of a s P atial P nase modulator and resolution elements 

phase contrast imaging is not disclosed in this document. of toe generated intensity pattern. 

A similar example of a phase contrast imaging system is According to the invention a method is provided for 
disclosed in "Array illuminator based on phase contrast", synthesizing an intensity pattern with low loss of e ectro- 
Applied Optics Vol. 27, No. 14, pp. 2915-2921 (1988). A ma S nehc en f uprising spatial modu ation of electro- 
method is disclosed of converting a wide beam of uniform 25 magneUc radiation with a spatial phase mask for modulation 
intensity into an array of bright spots without losses. Tlie °f thc °[. ™ ld * ni cfectromagnctic radiation by 

f *■ 1 u „ _ • i- „„.„., nf phasor values of individual resolution elements or the spatial 

input spatial phase mask constitutes a periodic array 01 Sf . 1 , • 1 t • . 

u j „ r - t u ..u„^ ~ • • nran „f f u„ phase mask, each phasor value being determined in such a 

phase dots with the phase value Jt, the remaining area or the r ' r & 

phase mask having the phase value 0. The phase modulation wa ^ a 

is converted into an amplitude modulation by Fourier trans- 30 1) the values of the Fourier transformed phasors attains 

forming lenses and a phase plate providing a phase shift of predetermined values for predetermined spatial 

the background signal by Jt. The method is limited to the frequencies, and 

implementation of periodic array configurations with the 2) the phasor value of a specific resolution element of the 

binary phase values 0 and jr. spatial phase mask corresponds to a distinct intensity 

It is well-known to use so-called "radiation focusators", 35 | evel of th e image of the resolution element in the 

i.e. computer generated holographic optical elements, for intensity pattern, 

spatial phase modulation of a light beam, e.g as disclosed in and a spatial phase filter for phase shifting of a part of the 

Special Issue on Computer Optics in the USSR, Optics and electromagnetic radiation, in combination with an imaging 

Lasers in Engineering, Vol. 15, no. 5 1991. However, such system for generation of the intensity pattern by interference 

elements are complicated to synthesize. Typically, they are 40 * toe image plane of the imaging system between the part 

synthesized in such a way that the desired image is formed of toe electromagnetic radiation that has been phase shifted 

in the Fresnel region or the Frauenhofer region. Thus, the b y | he phase filter and the remaining part of the electromag- 

intensity of a resolution element in the generated image is a netic radiation. 

function of several, typically all, phase values of the reso- Although, the present method is related to encoding of 

lution elements of the holographic optical element. 45 spatial phase masks in two spatial dimensions (planar 

Obviously, this complicates the design of a general purpose encoding), the principles of the method may be utilized for 

holographic optical element and advanced, very time con- P hase encoding in one to three spatial dimensions and/or in 

suming algorithms have to be applied. Further, the compli- th e temporal dimension. 

cated design of the holographic optical elements renders it The electromagnetic radiation may be of any frequency 

almost impossible to implement dynamically changeable 50 ran S e of tQ e electromagnetic spectrum, i.e. the gamma 

spatial phase modulators with such elements. frequency range, the ultraviolet range, the visible range, the 

tj • p ii 1 • I . e . , , . 1 1 infrared range, the far infrared range, the X-ray range, the 

It is a further disadvantage 01 holographic optical ele- . to ' it _ T ,„ ... . c * . J ~ ' 

.... c 0 ?/, * A t microwave ranee, the HF (high frequency) ranee, etc. The 

ments that a earner frequency is needed to separate dif- . , b . ' S , , . M • / 

c . , v t.r t jv u* u- • w • present method is also applicable to particle radiation, such 

fracted fight from non-diffr acted light resulting in an off-axis „„ v . . vv ,. " y 

system geometry and a need for a diffractive medium that 55 « electron radiat.on, neutron radiation, etc. 

can support these high frequency terms. Preferably, the electromagnet* ndi. .on » monochro- 

maUc or quasi-monochromatic so that the energy or the 
electromagnetic radiation is concentrated in a narrow fre- 
quency bandwidth. As the intensity pattern is synthesized by 

T# • . . . „ f » _ t • t„ „„ 60 interference of two electromagnetic waves emitted from a 

It is an object ot the present invention to provide an - . ^ *» . . , A it _ . 

apparatus of the above kind which apparatus is robust, c ° m ™° f™? of ^dromagnetic radiation but the phases 

compact, simple to design and relatively cheap to manufac- ? f whlch have been changed cufferenay, U is required that the 

ture r r ° j r frequency range of the emitted electromagnetic radiation is 

sufficiently narrow to ensure that the two waves of electro - 

It is another object of the present invention to provide an 65 magnetic radiation are coherent so that their superposition 

improved method and apparatus for phase contrast imaging generates the desired intensity pattern. If the frequency 

that take all terms of the Taylor's series: range is too broad, the two waves will be incoherent and the 
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phase information will be lost as superposition of non- 
coherent waves results in a summation of the intensities of 
the two waves. It is required that the difference between 
individual delays of electromagnetic radiation to be super- 
positioned is less than the wavelength of the radiation. This S 
is a relaxed requirement that allows the electromagnetic 
radiation to be relatively broad-banded. For example in the 
visible range a Xe-lamp or a Hg-lamp can be used as a light 
source in a system according to the present invention with 
the advantage compared to a laser light source that the 10 
speckle noise is reduced. The requirements of the spatial 
coherence of the electromagnetic radiation depend upon the 
space bandwith product of the corresponding system and 
how close the required system performance is to the theo- 
retically obtainable performance of the system. 15 

Preferably, the electromagnetic radiation is generated by 
a coherent source of electromagnetic radiation, such as a 
laser, a maser, a phase-locked laser diode array, etc. How- 
ever a high pressure arc lamp, such as a Hg lamp, a Xe lamp, 
etc, may also be used and even an incandescent lamp may 20 
be used as a source of electromagnetic radiation in a low 
performance system. 

A spatial phase mask is a component that changes the 
phase of an electromagnetic wave incident upon it The 
spatial phase mask may transmit or reflect the incident 25 
electromagnetic wave. Typically, the spatial phase mask is 
divided into a number of resolution elements each of which 
modulates the incident electromagnetic wave by changing 
its phase by a specific predetermined value. The predeter- 
mined values are assigned to each resolution element in 30 
different ways depending upon the technology applied in the 
component. For example in spatial light modulators, each 
resolution element may be addressed either optically or 
electrically. The electrical addressing technique resembles 
the addressing technique of solid-state memories in that each 35 
resolution element can be addressed through electronic 
circuitry to receive a control signal corresponding to the 
phase change to be generated by the addressed resolution 
element. The optical addressing technique addresses each 
resolution element by pointing a light beam on it, the 40 
intensity of the light beam corresponding to the phase 
change to be generated by the resolution element illuminated 
by the light beam. 

Spatial phase masks may be realized utilizing fixed phase 
masks, devices comprising liquid crystals and being based 45 
on liquid crystal display technology, dynamic mirror 
devices, digital micromirror arrays, deformable mirror 
devices, membrane spatial light modulators, laser diode 
arrays (integrated light source and phase modulator), smart 
pixel arrays, etc. 50 

A spatial phase filter is typically a fixed phase mask, such 
as an optically flat glass plate coated with a dielectric layer 
at specific positions of the glass plate. However, the spatial 
phase masks mentioned in the previous section may also be 
used for spatial phase filters. 55 

The imaging system maps the phase modulating resolu- 
tion elements of the spatial phase mask on the target surface 
of the synthesized intensity pattern. It may comprise a 
4f-lens configuration (two Fourier transforming lenses uti- 
lizing transmission of light or one Fourier transforming lens 60 
utilizing reflection of light) or a single imaging lens. 
However, any optical imaging system providing a filtering 
plane for the spatial phase filter may be applied in a phase 
contrast imaging system. 

In the method according to the present invention, the 65 
synthesized intensity pattern is generated by superposition 
of two electromagnetic waves in the image plane of the 



imaging system. The spatial phase mask changes the phase 
values of an electromagnetic wave incident upon it and the 
imaging system directs the electromagnetic wave with 
changed phases reflected from or transmitted through the 
spatial phase mask towards the spatial phase filter. The phase 
filter phase shifts a part of the electromagnetic radiation and 
the imaging system is adapted to superimpose in the image 
plane the phase shifted part of the electromagnetic radiation 
with the part of the electromagnetic radiation that is not 
phase shifted by the spatial phase filter. 

According to a preferred embodiment of the invention, 
the spatial phase mask is positioned at the front focal plane 
of a lens while the spatial phase filter is positioned in the 
back focal plane of the lens, whereby a first electromagnetic 
field at the phase mask is Fourier transformed by the lens 
into a second electromagnetic field at the phase filter. Thus, 
specific spatial frequencies of the first electromagnetic field 
will be transmitted through the spatial phase filter at specific 
positions of the phase filter. For instance, the energy of the 
electromagnetic radiation at zero frequency (DC) is trans- 
mitted through the phase filter at the intersecting point of the 
Fourier plane and the optical axis of the lens also denoted the 
zero-order diffraction region. 

It is presently preferred that the spatial phase filter is 
adapted to phase shift the DC-part of the electromagnetic 
radiation and to leave the remaining part of the electromag- 
netic radiation unchanged or, alternatively, to leave the 
DC-part of the electromagnetic radiation unchanged and to 
phase shift the remaining part of the electromagnetic radia- 
tion. The last alternative is preferred when the energy level 
of the DC-part of the electromagnetic radiation is so high 
that the phase shifting part of the phase filter will be 
destroyed by it. For example in laser cutting, the DC level 
of the laser beam can be so high that a phase shifting dot 
positioned at the intersecting point of the DC part of the laser 
beam at the phase filter would evaporate. It is also possible 
to block the electromagnetic radiation (no transmittance) in 
the zero-order diffraction region, however, the DC energy of 
the radiation is then lost. 

Below, an expression of the intensity of the synthesized 
intensity pattern as a function of the phasor values <Kx,y) of 
the phase mask, when the DC-part of the electromagnetic 
radiation is phase shifted, is deduced. 

Electromagnetic radiation incident on the spatial phase 
mask can be described by a function A(x,y), where A(x,y) is 
a complex number (amplitude and phase) of the incident 
field on the point (x,y) of the spatial phase mask. At the point 
(x,y), the spatial phase mask modulates the phase of the 
incident radiation with a value <(>(x,y) so that the field after 
reflection by or transmission through the spatial phase mask 
may be described by the function A(x,y)*e ,vt>c ^>, e*^ 
being the phasor value of the point (x,y) of the spatial phase 
mask. As A(x,y) preferably is a constant value over the entire 
surface of the spatial phase mask, the term is left out of the 
following equations for simplicity. 

The expression of the electromagnetic radiation incident 
on the spatial phase filter may now be separated into an 
AC-term and a DC-term. If the DC-term of the field is 
denoted a, the AC-term of the field is given by the term 
e^c^yj-a. As the spatial phase filter changes the phase of the 
DC-part of the electromagnetic radiation by 6, the intensity 
of the synthesized intensity pattern at the image plane of the 
imaging system is given by: 

/(x^O-k^^+a^-l)! 2 (7) 

wherein (x',y') is the coordinates of the image of the point 
(x,y) of the spatial phase mask formed by the imaging 
system in the image plane. 
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It should be noted that the second term of the equation is 
a complex number that adds to the phasors e t<1,(x,y) of the 
spatial phase mask and may be interpreted as a contrast 
control parameter for the synthesized intensity pattern I(x', 

y). 

According to a preferred embodiment of the invention, 
the average value of the phasors is adjusted in order to 
control the range of intensity levels. 

Instead of phase shifting the DC-part of the electromag- 
netic radiation, it is also possible to synthesize a prescribed 
intensity pattern by phase shifting other parts of the elec- 
tromagnetic radiation by adapting the spatial phase filter to 
phase shift electromagnetic radiation incident upon one or 
more arbitrary regions of the phase filter and leaving the 
phase of the remaining part of the electromagnetic radiation 
unchanged and then superimposing the two parts of the 
electromagnetic radiation. The corresponding mathematics 
and the corresponding design procedures for the spatial 
phase mask and spatial phase filter will of course be more 
complicated than for the method described in the previous 
section. 

A simple example of phase shifting a part of the electro- 
magnetic radiation of a spatial frequency different from the 
zero frequency is provided by moving the DC-part of the 
electromagnetic radiation to another spatial frequency in the 
Fourier plane (identical to the plane of the spatial phase 
filter) utilizing an optical component with an appropriate 
carrier frequency (i.e. a grating or a prism) or, preferably, 
encoding the function of a grating or a prism into the spatial 
phase mask, and adapting the spatial phase filter to change 
the phase of the electromagnetic radiation at this spatial 
frequency and to leave the phase of the remaining part of the 
electromagnetic radiation unchanged. 

According to another preferred embodiment of the 
invention, the phase mask is not positioned in the back focal 
plane of the lens but in the Fresnel region of the lens instead. 
In this case, the electromagnetic field at the phase filter will 
be given by a Fresnel transformation of the electromagnetic 
field at the spatial phase mask. This further complicates the 
mathematics and the design procedures, for example the 
term a in equation (7) has to be substituted by the value of 
the Fresnel transformation at the point(s) of phase changes 
of the phase filter. However, the Fresnel transformation may 
be calculated from a Fourier transformation by multiplica- 
tion of the phasor values of the spatial phase mask by a 
quadratic phase factor followed by a Fourier transformation. 

It is an important aspect of the present invention that each 
intensity level of the synthesized intensity pattern for each 
resolution element may be generated by at least two different 
phasor values of a resolution element of the spatial phase 
mask. 

For example, when the spatial phase filter phase shifts the 
DC-part of the electromagnetic radiation, it will be shown 
later that, advantageously, the average a of the phasors of 
the resolution elements of the phase mask should be equal to 
l A and the value of the phase shift 0 should be equal to Jt. In 
this case, the intensity of the synthesized image pattern at the 
image (x\y') of the resolution element (x,y) will be given by: 

It is seen that complex conjugate phasors (values of <|> of 
opposite sign) result in identical intensity levels I(x',y'), It 
can be shown that for any value of the modulus of the 
average of the phasors |a|, two phasors exist that will 
generate identical intensity levels of the synthesized inten- 
sity pattern. 

Further, if the spatial phase filter phase shifts parts of the 
electromagnetic radiation different from the DC-part, the 
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phasor value that generates a specific intensity level will 
depend on the position of the resolution element in question, 
i.e. the phasor value and the position of the resolution 
element with that phasor value together define the intensity 

5 level at the image of the resolution element in the synthe- 
sized intensity pattern. Still, it is true that for each resolution 
element of the spatial phase mask, each intensity level of the 
synthesized intensity pattern may be represented by one of 
two different phasors of complementary phase values. 

10 This freedom of being able to select, for each intensity 
level to be generated and for each resolution element of the 
spatial phase mask, one of two phasors is used to control the 
phase of the Fourier transform of the phasors at specific 
spatial frequencies by selection of phasors with appropriate 

15 phase values to ensure two intervals of biunique functional 
dependence between phasor values and corresponding inten- 
sity values. 

This freedom of choice of phasors may be utilized to 
select phasors of neighbouring resolution elements of the 

20 spatial phase mask with a maximum difference between 
them, thereby generating an electromagnetic radiation emit- 
ted from the phase mask with a maximum content of high 
spatial frequencies which will generate a good separation of 
the DC part of the electromagnetic radiation from its AC 

25 part. However, any other strategy of selecting between two 
possible phasor values of each resolution element may be 
chosen to generate a desired spatial frequency content of the 
electromagnetic radiation. 
Preferably, the phase of the Fourier transform of the 

30 phasors at specific spatial frequencies is adjusted in order to 
control whether the relation between each phasor and the 
corresponding intensity level is a mono tonic increasing or a 
mono tonic decreasing function. 

Below, a set of different methods are described that are 

35 provided according to the present invention for adjustment 
of the modulus of the Fourier transform of the phasors at 
specific spatial frequencies to attain a prescribed value. If 
convenient, the methods may be combined. 
According to one of the methods, the individual phasors 

40 of the resolution elements of the phase mask are adjusted by 
a constant value until the desired value of the modulus of the 
Fourier transform of the phasors at specific spatial frequen- 
cies is attained while maintaining prescribed relative inten- 
sity levels between intensities of resolution elements of the 

45 intensity pattern, i.e. iteratively. 

According to another method, the individual phasors of 
the resolution elements of the phase mask are adjusted 
utilizing histogram techniques known from image process- 
ing. A histogram is a bar chart showing the number of 

50 resolution elements of the synthesized intensity pattern with 
a specific intensity value as a function of the intensity value. 
Any histogram technique, such as histogram equalization, 
adapting the histogram to a predetermined distribution, etc., 
may be used iteratively until the modulus of the Fourier 

55 transform of the phasors at specific spatial frequencies attain 
the prescribed value. 

According to yet another method, the phasor pattern of the 
phase mask is spatially scaled in order to adjust the modulus 
of the Fourier transform of the phasors at specific spatial 

60 frequencies. 

According to still another method, the modulus of the 
Fourier transform of the phasors at specific spatial frequen- 
cies is adjusted utilizing half tone coding techniques, such as 
raster techniques, area ratio modulation, spot diameter 

65 modulation, etc. 

It is seen from the description above that the intensity 
levels may differ from one synthesized intensity pattern to 
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the next as a consequence of the adjustments of the modulus expander (3) into a plane light wave of uniform intensity and 
of the Fourier transform of the phasors at specific spatial directs it towards a spatial phase mask (4). The light beam 
frequencies. Thus, it is preferred to control the power of the is transmitted through the spatial phase mask (4) and a 
radiation source in dependence of the intensity range of the Fourier transforming lens (5). The spatial phase mask is 
intensity pattern so that a sequence of different intensity 5 positioned in the front focal plane of the lens (5) and a spatial 
patterns show uniform intensity levels. P hase fllter W * positioned in the back focal plane of the 

According to a preferred embodiment of the invention, lens ( 5 ) that 15 also the front focal plane of a lens (7). The 
the shape of the phase filter is adapted to match the spatial Courier transforming lenses (5, 7) need not have identical 
frequency content of the phasors of the spatial phase mask, focal knglhs. Mkrent focal lengths lead to a magnification 
eg. to optimize the desired separation of the part of the 10 ratio ^ r ™} j* om one - ^ ^phase mter (6) phase shifts the 
electromagnetic radiation to be phase filtered from the zero order diflfrachon part (8) of the light phase modulated 
remaining part of the electromagnetic radiation. b V the s P atial P hase W- ™ e ^^f 2 ?* mteQS "y 

It is within the scope of the present invention that the pattern is generated in the back focal plane (9) of the lens (7) 
imaging system further comprises zooming means for vari- and a dynamic focusing system (10) images the synthesized 
able scaling of the synthesized intensity pattern. The zoom- 15 intensit y P attcrn onto a focusin S P lanc ( n >* 
ing of the imaging system may be dynamically controllable, The optical system is controlled by a computer (12). The 
e.g. in response to the scaling of the pattern of phasor values computer (12) comprises interface means for addressing 
of the phase mask. eacn °f me resolution elements of the phase filter (4) and 

According to the present invention, the power of the transmitting a phasor value to the addressed resolution 
radiation source may be controllable in response to the 20 element. Further, the computer (12) comprises laser control 
spatial scaling of the pattern in the phase mask and/or the means for controlling the power of the laser (2) and imaging 
zooming of the focusing system. control means for controlling the focusing and the image 

In order to provide a compact and integrated system ratio of the dynamic focusing system (10). The computer 
according to the present invention, the optical function of a (12) also comprises input means, such as a keyboard, a 
Fourier-transforming lens is encoded into the phasors of the 25 diskette drive, an optical disc drive, a network interface, a 
spatial phase mask. The Fourier transforming lens may be modem, etc, for receiving an image pattern to be synthesized 
refractively or diffractively encoded into the phase mask. by the system (1). From the received image pattern, the 

Similarly, the optical function of an output lens may be computer is adapted to calculate phasor values to be trans- 
encoded into the phase filter either refractively or diffrac- mitted to the resolution elements of the phase mask, e.g. 
t * ve jy 3 q based on a histogram technique as described herein. 

Further, a compensation may be encoded into the phasor Optionally, the phase shift of the phase filter (6) is adjustable 
values of the spatial phase mask so that part of the electro- ^ controllable by optional phase control means of the 
magnetic radiation modulated by the phase mask has a computer (12) which may be further adapted to adjust the 
substantially flat intensity profile in the image plane. With- P^se shift, e.g. utilizing equation (18). 
out this compensation, part of the electromagnetic radiation 35 FIG. 2 shows a 2f phase contrast imaging system (20). A 
modulated by the phase mask will have a flat profile with laser (21) emits a light beam which is expanded by a beam 
perturbations resulting from the phase filtering superposi- expander (22) into a plane light wave of uniform intensity 
tioned upon it. This may cause "ringings" (oscillations) at and directs it towards a spatial phase mask (23) and a 
the edges of the synthesized intensity pattern. polarization beam splitter (24) and a quarter-wave plate (25). 

According to another preferred embodiment of the 40 The polarization beam splitter (24) and the quarter-wave 
invention, the source of electromagnetic radiation comprises plate (25) allows beam-splitting of fight of a specific linear 
one or more light sources of different wavelengths corre- polarization without the power loss associated with conven- 
sponding to three different colours, such as red, green and tional beam-splitters due to splitting of the beam in both 
blue, for generation of intensity patterns of arbitrary colours. directions of transmission through the beam-splitter. After 
Further, several independent systems each one illuminated 45 transmission through the polarization beam splitter (24) and 
by its own wavelength can be combined into a single the quarter- wave plate (25), the light beam is transmitted 
multi- wavelength system. through a Fourier transforming lens (26) and is reflected 

from a spatial phase filter (27). The spatial phase mask (23) 
BRIEF DESCRIPTION OF THE DRAWINGS is positioned in the front focal plane of the lens (26) and the 

HG. 1 shows a 4f optical system for phase contrast 50 spatial phase filter (27) is positioned lin the back focal plane 

imagine of the lens ( 26 )* The pnase filter ^ phase shlfts 1 zero 

JL ' , v , , , , t „ t order diffraction part (28) of the fight that is phase encoded 

FIG. 2 shows a 2f optical system for phase contrast * v \ m . * , r . , . . 

r 3 r by the spatial phase mask (23). The synthesized intensity 

imaging, pattern is generated in the back focal plane (29) of the lens 

FIG. 3 shows a if optical system for phase contrast 55 (2fi) afld ft dynamic focusing systcm (30) images thc syn . 

imaging, thesized intensity pattern onto a focusing plane (31). As 

FIG. 4 shows (A) off-axis read-out of reflective SLM and described for the system shown in FIG. 1, the system (20) is 

(B) on-axis read-out of reflective SLM. controlled by a computer (32). 

FIG. 5 shows schematically an example of a prescribed fig. 3 shows a if phase contrast imaging system (40). A 

intensity pattern in ID. 60 i aser (41) emits a light beam which is expanded by a beam 

FIG. 6 shows schematically the resulting phase encoding expander (42) into a plane light wave of uniform intensity 

corresponding to FIG. 5. and directs it towards a spatial phase mask (43). The light 

beam is transmitted through the spatial phase mask (43) and 

DETAILED DESCRIPTION OF PREFERRED an ^ forming leDS (44) . A phase filter (45) positioned in 

EMBODIMENTS 65 the back focal pUne of the lens ^ phase &Mis the ^ 

FIG. 1 shows a 4f phase contrast imaging system (1). A order diffraction part of the light phase encoded by the 
laser (2) emits a light beam which is expanded by a beam spatial phase mask (43). The synthesized intensity pattern is 
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generated in the image plane (46) of the lens (44) and a 
dynamic focusing system (47) images the synthesized inten- 
sity pattern onto a focusing plane (48). As described for the 
system shown in FIG. 1, the system (40) is controlled by a 
computer (49), 

FIG. 4 shows details of (A) an off- axis read -out of a 
reflective phase mask (50) (or a spatial light modulator) and 
of (B) an on-axis read-out of a reflective phase mask (51) 
with a beam splitter (52). Both configurations (A, B) may be 
utilized in the systems shown in FIGS. 1-3. 
Phase Encoding For DC Phase Filtering 

In the following an example of encoding a spatial phase 
mask and a spatial phase filter will be given based on a 
system filtering in the DC-frequency range. The system 
chosen in this example is based on a 4-f lens configuration 
as shown in FIG. 1 and illuminated by electromagnetic 
radiation in the visible frequency domain, hereafter simply 
denoted as light radiation. 

Assuming that the Uluminating light is monochromatic 
and has a substantially flat amplitude profile we obtain the 
following spatial amplitude distribution emitted from the 
spatial phase mask: 



a{x, y) = rect (£^ ^jexp(*0(;t, y)) 



A* Ay//' 



exp(i<p(x, y))dxdy 



(10) 



Subsequently the AC-term is expressed by: 

Aa(x, y) = exp<t0(jr, y))- J J exp(i^(jc, y))dxdy 



(11) 



The separation of ct(x,y) into a spatially invariant 
DC-term and a spatially varying AC-term is an important 
point and will be used throughout the remaining part of this 
example, especially in the description of the spatial filtering 
procedure. 

The spatial filter utilized in this example is chosen as a 
circular phase contrast filter (different transverse shapes can 
also be used) centered around origo in the the spatial 
frequency domain, denoted by coordinates (f*,f y ): 



T(/ r )=l+(exp(*0)-l)cir< 



(12) 



where f ^ y% 2 +f y 2 denotes radial spatial frequency and Af^ 
describes the size of the circular (circ) phase filter. 

In the spatial frequency domain (the filtering plane) the 
Fourier transformation (5) of the spatially modulated light 
radiation from the spatial phase mask is present. The filter- 
ing operation on the Fourier transformed light radiation 
performed by the spatial phase contrast filter can be 
expressed as a simple point-by-point multiplication proce- 



12 



dure. Subsequently the spatially filtered light is inverse 
Fourier transformed (S _1 ) by the second Fourier lens 
(Fourier transformation and reflected output coordinates) 
and the resulting spatial amplitude distribution in the image 
plane (with coordinates (x',y')) can accordingly be written 
as: 



15 



20 



o(*\ /) = a(*\ /) + (exp(ifi) - l)T- l \r{a{x, y))cir^ A jj 
a^,/) + ZF(«p(W)-l)rect^, 
= [expWt*', /)) + <F( C xp(/0)- Direct^, 



(13) 



Within the illumination-icgion, (x\y')e9?, outlined by 



(9) 



one obtains: 



where a(x,y) =exp(i(()(x,y)) represent the spatially encoded 
phasor values and AxAy is the area of the input phase 
modulating spatial light modulator. 

It turns out to be convenient to separate a(x,y)_into two 
terms describing a spatially invariant DC- value, a, and a 
spatially varying AC-contribution Aa(x,y). The DC-value 
can be found as: 



30 



C*\ St s 1 + 4|*|sjnQ[|c?| S iii^] - sin^ y')+ j] 



Requiring that |o(x 0 ',y 0 ')| 2 =0 corresponding to complete 
darkness as the lowest intensity level in regions (x 0 \yJ)ESiJ 
implies: 



35 



1 + 



(15) 



where the abbreviation ^oH^O^y©*) nas Deen used. 
The solutions to Eq. (15) are given by: 



W = 



45 



sin^o -fa + ± ^sin 2 ^ -fa + - 



(16) 



50 



The requirement o<|a]<l implies that: 
0 = n - 2^ a - fa) + plr, /> = 0, ±1, ±2, .. 



(17) 



55 leading to 



±1 



2sin(-] 



(18) 



60 



where the +sign is for 8-values in the interval: 



6 5 e \y 5 i\^^ 



(19) 
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and the -sign is for 6-values: 



9e \y~3\ +Podd2jr 



(20) 



The corresponding interval for (fe-<t> 0 ) i g: 



7i n 
3 ; ~ 3 



(21) 



Inserting the expression for |a|, one obtains the simple 
intensity expression: 



/)l 2 = 2|i T sin(fc -*(*\ /)+■ jj] 



(22) 



where 



JJexp(i<p(x, y))dxdy^ Lxky\a\cxv{ty m ) 



(23) 



25 



The phase-only transformations imply that energy is 
conserved: 



J J |o(*\ /fdx'dy' = J JW ytfdxdy = *xA.y 



(24) 



A special case: 

The most convenient choice for a is: a=Vz (implying that 
e«=n;+P eve . M 2jt), so that the output intensity can be described 
as: 



M*;y0r=2[i-cos(<K*',>''))] 



In this case the phaseointensity mapping is described by 
the intervals [O;*]— [0;4]. 

By setting a.= x A one obtains the following requirement to 
the phase function ((>(x,y): 



J Jcos(<p(x t y))dxdy- 
J J s'm(<(>{x, y))dxdy = 0 



AxAy 



(26) 



Inserting the expression for ^(x'^y')! 2 in Eq. (24) yields: 

2jJ[i- cos^y, y'))]dx'dy J = AxAy (2?) 



in accordance with the first of the integral expressions in Eq. 
(26). 

Encoding procedure: 

A given intensity distribution (image) ^(x^y 1 )! 2 is desired at 

the output side of the optical setup. 
Pixellation of the image, that is generally represented in the 

greyscale range: [0;gmax], provides the relation: 



30 



35 



(25) 



45 



55 



14 



The histogram for the desired image |o(ij)| 2 is adjusted (adj) 
within the greyscale range [0;gmax], so that the previous 
point is fulfilled: 

10 |o(ij)| 2 -|o(ij)U/. 

The phase values can now be calculated as: 



${i t j) = arccos 



y gmax ) 



As before pixellation provides the relation: 

20 ^sm(0(i,y)) = O. 



50 



The previous point can now be fulfilled by complex conju- 
gating half the input pixels having the same phase value 
in the phase histogram. 

The phase conjugate phase flipping provides a valuable tool 
(an extra degree of freedom) for manipulating the spatial 
frequency content in order to optimize the separation of 
low and high frequency terms at the filter plane. 

The scheme is robust to constant phase errors across the 
input spatial phase modulator, since Eq. (22) is a function 
of the difference: <|>£-<Kij), on ly- Furthermore, small 
variations in the individual pixel phase values do not 
introduce any detrimental effects because the average 
value a, is a result of a very large phasor sum. 

If the desired intensity distribution is too small to include all 
energy, that is, the histogram is scaled to maximum and 
the left hand side of Eq. (24) is still smaller than the right 
hand side, then the input phase object can be scaled until 
Eq. (24) is fulfilled. In order to obtain a scale invariant 
output intensity level a dynamic focusing system is 
needed. Similarly, intensity invariance can be obtained by 
controlling the radiated power from the light source. 
Alternatively, one can ignore the residual background 
illumination and obtain intensity levels with a gain factor 
of. 9~ (background constant equal to 1") for narrow 
generally shaped line structures (e.g. Eq. (14)). 

EXAMPLE 1 

A very simple example illustrating the individual steps in 
the above procedure will be given below. To simplify the 
example it will be considered in one dimension only. The 
starting point for encoding the spatial phase mask in this 
example is based on the following parameters: 



60 



w = 0.5 



(28) 



# pix {t 
gmax = 4 



: 14 



Consider the pixellated 3-step function shown in FIG. 5 to 
be synthesized in the image plane as an intensity distribu- 
65 tion. From the above choices of parameters one obtains the 
simple relation between phase values in the spatial phase 
mask and the image intensity values: 
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K0| 2 =2[l-cos(<t>(0)] (29) 

To proceed from here it necessary to calculate the accu- = - (4exp(i0) + 2exp(i0.97) + 

mulated intensity 3 « P (il.4S) + 2ex P (-i0.97) + 3ex P (-il.45)) = l /2 

5 

£>(0I 2 

' General Phase Correction Procedure Integrated With the 

Phase Encoding 

in the image to be synthesized. The accumulated intensity is [ n £q. (14) we obtained an analytic relation between the 

easily calculated from an image histogram where the x-axis phase va i ues m the spatial phase mask and the resulting 

represents greylevel value and the y-axis represents the . • , , 

. • i • *u ■ * * ill tj., intensity distribution, within the region (x',y')e9l ': 

amount of pixels in the image at a given greylevel value. By J * ^ v " ' 

use of a histogram 



£>(i)l 2 



(36) 



The analysis leading to the above relation was based on 

is simply found as the weighted sum of all greylevel values the assumption that |o| is a constant value within the 

(x-axis) multiplied by their pixel counting (y-axis). This 20 91' -domain. In other words, the following approximation was 

describes, so to speak, the "weight" of the image. In this applied: 
simple example histogram calculations are not needed since 

we only have 3 greylevels with well-defined separations. , ( f r \\ ( * y \ (37) 

The value for the accumulated intensity has to obey the r^rw*. >)^xjrJJ affn ^£» 



equality: 



max-?/* (32) 



£ _ 8max n P ix {iu() = # pix^ = 14 (30) However, for certain spatial filter parameters the lefthand 

/ 4 side of this expression will not be a space invariant constant 

30 value throughout the whole St'-domain but will instead 
From FIG. 5 we obtain: manifest slowly variations/oscillations. This will introduce 

small errors in the final superposition between the phase 
2 , , • , o on filtered DC-value and the direct propagated AC-signal. In 

o(/)r = 4 pixels -0 + 4 pixels -(0.5max) + 6 pjxcls-mav= %max \ Jlf . , . , . . * , . t 

. order to circumvent this problem a technique is needed that 

35 can counteract the distortions by use of phase -only encoding 
in the components already present in the system. In what 
So that the value for max can be estimated to be: 7 follows a procedure for integrating predistortion that coun- 

teracts the above mentioned distortions will be described 
that is purely based on modifying the phasor values in the 

Hie corresponding adjusted intensity levels, |o(i) * are *" s P atial P hase mask at the in P* side of f * e s y stem ; 

therefore: %, 7 /s and 0, These values can now be utilized to metbod can also counteract other types of distortions mher- 

calculate the phase values of the spatial phase mask from the ent in a practical implementation of the system, 

relation* Furthermore, the method can be applied in systems filtering 

45 at other spatial frequencies than DC. 

( 2\o(i)&j )_ ( MOI^) 03) Procedure: 

m ~ " ~^T) " ^y ' 2 J When encoding the input phase function it is helpful to 

have a reverse equation, expressing the input phase distri- 
bution as a function of an adjusted (electronic) image 

where from we obtain the three phase values: 1 .45 rad. 0.97 50 grey . level distribution, \ slm , addressing the input spatial light 

rad. and 0 rad. modulator: 

The last step needed in order to encode the spatial phase 

mask is that the following equality is fulfilled: 4 , . r g^r ns) 

* 1 +4|ffU\ /)|sii^-||<F(*\ /)|sin^- J - 
y>(*C0) = 0 (34) " / 0X\ 

^<My,/)-0(x\ /) + -]] 



Since we have the choice to use complex conjugate 
phasor values (two phasors giving the same intensity level) 
many approaches can be taken from here. A simple approach 
is to flip every second phasor with its complex conjugate 
value as shown in FIG. 6. The final phase values used in the 
phase masks are accordingly: ±1.45 rad. ±0.97 rad. and 0 
rad. 

As the last step we can check whether the criteria: a«H, 
is actually fulfilled with the chosen phasor encoding: 



where it has been taken into account that a(x',y^ is not 
60 considered as a constant but manifests a smooth oscillating 

behaviour within the optical image domain. The maximum 

value of \ slm is denoted gmax. 

Now, one can derive a formula for the grey-level 
65 correction, AI j/m (x , > y'), that one needs to apply in order to 

encode a phase function that compensates for the spatial 

variation of the average phase value a(x , ,y f ): 
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( /) , M . ^ . (0) (39) 

a 1 + 4|a(y, /)|sud - 

grnax \2) 

<M*\ y ) = acos 1 

V grnax J 



J(*. y) = 



/(*. y) 



(46) 



where the second relation has been derived from the first by 
setting a=VS and 9=jt. 

By inserting the second relation in the first expression one 
gets: 



f — - 

V2IW. 



(40) 



grnax 



This formula is however not directly useful because it is 
related to the histogram adjusted grey-level distribution 
denoted by \ slm . 

One needs a formula that relates the above correction term 
to the original input grey-level distribution I(x,y) that has not 
been modified by histogram adjustments. This is important 
since the effect of the grey-level corrections also have to be 
incorporated in the procedure of histogram adjustments. 

The histogram scaling gives: 



/(*. y) = 



htm. IX 



- W*> y) 



(41) 



where \ max and \ slmjnax 



are the maximum grey-level values 
occurring in the original and the adjusted electronic grey- 
level distributions respectively. 

Similarly, one can apply this relation to the intensity 
correction term M slm and obtain: 



7U, y) = /<*, y) + A/(* y) = ( W*, y) + AW*, y)) m 



.ID 



resulting in: 

2\cS{x, y)\\ Wm«v 2 J J 



In order to have enough dynamic range in grey -levels for 
the correction term one can derive an inequality from the 
above relation by using the fact that i max = grnax: 



215^1 V ishn.nwV 2} J 



215^111 1 -gmax 



Since the first term is the dominating term in the expres- 
sion for the intensity correction it will in practice be suffi- 
cient just to have the much simpler corrections: 



Proposed Applications: 

Laser machining, marking, branding, trimming, hardening, 
scribing, labeling, welding and cutting on two- and three- 
io dimensional surfaces especially by use of C0 2 and 
Nd:YAG laser based systems. The main advantage is that 
energy is not absorbed in the system (thereby preventing 
damage of the optical hardware) and this non absorbed 
energy is instead utilized to increase the intensity level of 
15 the desired light distribution in the image plane. High 
power can be delivered to selected regions on a work 

(1^2 piece simultanously. 

|5f(V, /)| - -J Efficient and dynamic spot-array generators based on phase 

contrast imaging. In order to provide bias or holding 
20 beams for arrays of optoelectronic elements, such as 
bistable elements, photonic switches and smart pixels. 
Generation of structured light (lossless) for machine vision 
applications. E.g. periodic and skew periodic mesh grid 
illumination that can be updated in parallel. 
25 Photolithographic applications (laser 3D direct writing in 
parallel without the need for sequential scanning). E.g. 
high power laser direct writing of waveguides in 
Ge -doped silica. 
Spatial light intensity modulation in general by use of pure 
30 phase modulation (radiation focusators). 
Laser beam shapinq (dynamic). 

Highly efficient parallel image projection without the need 

for a laser scanning device. 
Dynamic Infrared Scene Projection (DIRSP). 
35 Exposure device for grating and mask production. 
UDAR applications. 
Laserprinting in parallel. 
Lasershow applications. 
Atmosphere research. 
40 I claim: 

1. A phase contrast imaging method of synthesizing an 
intensity pattern I(x',y') of an image, comprising the steps of 
pixellating the intensity pattern I(x\y') in accordance with 
the disposition of resolution elements (x,y) of a spatial 
45 phase mask (4, 23, 43) having 

a plurality of individual resolution elements (x,y), each 
resolution element (x,y) modulating the phase of 
electromagnetic radiation incident upon it with a 
predetermined phasor value e*****^, 
radiating electromagnetic radiation towards the spatial 

phase mask (4, 23, 43), 
Fourier or Fresnel transforming the modulated electro- 
magnetic radiation, 
phase shifting in a region of spatial frequencies compris- 
ing DC in the Fourier or Fresnel plane, the modulated 
electromagnetic radiation by a predetermined phase 
shift value 6 in relation to the remaining part of the 
electromagnetic radiation, and 
forming the intensity pattern by Fourier or Fresnel 
transforming, respectively, the phase shifted Fourier or 
Fresnel transformed modulated electromagnetic 
radiation, whereby each resolution element (x,y) of the 
phase mask (4, 23, 43) is imaged on a corresponding 
65 resolution element (x',y') of the image, 

calculating the phasor values e 1 *^ of the phase mask (4, 
23, 43) and the phase shift value 0 in accordance with 



1\ 2 \ (43) 



55 



(44) 



(45) 60 
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/(*'.y>k*"- y) +a(e B -i)P b) adjusting the individual phasors e**^ of the resolution 

for selected phase shift values 9, a being the average of the S^uef ^ ***** ^ **' 43) by hist0sram 

phasors e'* (:c,y) of the resolution elements of the phase mask . A . ;} ' x . , / A( v^ A r , 

AX\ c ) spatially scaling the phasor e VK pattern of the phase 

K% ~> } > , , . , r _ s mask (4, 23, 43), and 

selecting, for each resolution element, one of two phasor d) ^ ^ toQe ^ techniques. 

values which represent a particular grey level, and 9 A method according to claim x> comprising the 

supplying the selected phasor values e***^ to the reso- ste p of controlling power of the electromagnetic radiation in 

lution elements (x,y) of the spatial phase mask (4, 23, response to the intensity range of the intensity pattern. 

43). 1Q 10. A method according to claim 1, wherein each phasor 

2, A method according to claim 1, wherein the step of e***'^ of the phase mask (4, 23, 43) is selected from a set of 

calculating the phasor values comprises two determined phasors with complementary phasor values 

setting the synthesized intensity of at least one resolution e* 1(jc ' y) and e**' 2 ^ in such a way that a specific spatial 

element (x 0 ',y 0 ') of the intensity pattern to zero, and frequency distribution of the intensity of the electromagnetic 

calculating the phasor values e^> of the phase mask (4, a5 ra ^ n in ! he / 0UT ™ or Fres ? el Plane is attained. 

23 43) in accordance with method according to claim 10, wherein the phase 

' ' <K x >y) °f phasors e 1 *^ of adjacent resolution elements 

alternates between the two possible complementary phasor 
values e* 1 **^ and e'* 2 ^. 



1*1 * 



//•■ 

AxAy 



2 | s j n i| 12. A method according to claim 10, wherein the phasors 

I 2 I 20 e** 1 ^ and e l * 2(;Ciy) are complex conjugated. 

r , q^-i 13. A method according to claim 1, further comprising the 

i{x\ /) = 2 1 TsiiJ fa - fix 1 , /) + - jl s tep of phase shifting at selected spatial frequencies consti- 
tuting a region that is shaped to match the spatial frequency 

_ content of the phasors e**^ of the spatial phase mask (4, 

for selected phase shift values 0, being the phase of a. 25 23, 43). 

3. A method according to claim 2, further comprising the 14. A method according to claim 1, wherein modulus of 
step of selecting the phase shift 0=Jt, selecting |a|= 1 /4, and the average value |ct| of the phasors e 1 ^'^ ranges from 0.1 
calculating the phasor values e**^ of the phase mask (4, to 0.9. 

23, 43) in accordance with 15. A method according to claim 14, wherein modulus of 

30 the average value |a| of the phasors c^ x ^ ranges from 0.25 

/C*\ /) = 2[l - costly, /))] 10 °* 75 - 

16. A method according to claim 14, wherein modulus of 

sin(^(x, y))dxdy = 0. the average value |a| of the phasors e'*^ ranges from 0.4 

to 0.6. 

35 17. A method according to claim 14, wherein modulus of 
the average value |a| of the phasors e l0C *' y) is approximately 
0.5. 

4. A method according to claim 1, further comprising the lg A me thod according to claim 1 wherein the phase shift 
ste P s of 6 ranges from Jt/4 to 7n/4. 

moving the DC-part of the electromagnetic radiation to a 4 o 19. A method according to claim 1, wherein the phase 

second part of the Fourier or Fresnel plane, and shift 6 ranges from Jt/2 to 7>txJ2. 

phase shifting the Fourier or Fresnel transformed modu- 20. A method according to claim 1, wherein the phase 

lated electromagnetic radiation at the second part of the shift 0 ranges from 3tt/4 to 5?i/4. 

Fourier or Fresnel plane by 6 in relation to the remain- 21. A method according to claim 1, wherein the phase 

ing part of the electromagnetic radiation. 45 shift 0 is approximately ji. 

5. A method according to claim 4, wherein the step of 22. Amethod according to claim 1, further comprising the 
moving the DC-part of the electromagnetic radiation com- step of zooming the image for scaling of the intensity 
prises utilization of an optical component, such as a grating, pattern. 

a prism, etc, with an appropriate carrier frequency. 23. Amethod according to claim 22, wherein zooming of 

6. A method according to claim 4, wherein the step of so the image is dynamically controllable. 

moving the DC-part of the electromagnetic radiation com- 24. Amethod according to claim 22, wherein zooming of 

prises encoding the function of an optical component, such the image is controllable in dependence of the scaling of the 

as a grating, a prism, etc, with an appropriate carrier phase mask (4, 23, 43). 

frequency, into the spatial phase mask (4, 23, 43). 25. A method according to claim 22, further comprising 

7. Amethod according to claim 1, further comprising the 55 the step of controlling power of the electromagnetic radia- 
step of adjusting the modulus of the Fourier transform of the tion in response to the spatial scaling of the pattern in the 
phasors e'*^ at specific spatial frequencies in order to phase mask (4, 23, 43) and/or the zooming of the image, 
control the range of intensity levels of the synthesized 26. A method according to claim 1, wherein the step of 
intensity pattern. phase shifting comprises utilization of a spatial light modu- 

8. A method according to claim 7, wherein the step of 60 lator. 

adjusting the modulus of the Fourier transform of the 27. Amethod according to claim 1, further comprising the 
phasors e 1 *^ at specific spatial frequencies comprises at step of encoding the optical function of a Fourier- 
least one of the following measures: transforming lens into the phasors e***'^ of the phase mask 
a) adjusting the individual phasors e'*^ of the resolution (4, 23, 43). 
elements of the phase mask (4, 23, 43) maintaining 65 28. Amethod according to claim 1, further comprising the 
prescribed relative intensity levels between intensities step of encoding the optical function of an output lens into 
of resolution elements of the intensity pattern, the phase filter (6, 27, 45). 
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ing 



20 



e**^ and being interlaced with the remaining rows or 
columns having the pbasor value e* 3 . 
31. A system (1) according to claim 30, further compris- 



means for pixellation of the intensity pattern I(x' 7 y l ) in 60 

accordance with the elements (x,y) of the spatial phase 

mask (4, 23, 43), 
means for calculating the phasor values e^^ of the 

phase mask (4, 23, 43) and the phase shift value 9 in 

accordance with equation 

/(^'^O-k^'^C^-i)! 2 , 



65 



22 



10 



15 



29. A method according to claim 1, wherein the step of 
radiating electromagnetic radiation comprises radiation of 
electromagnetic radiation of different wavelengths corre- 
sponding to three different colours, such as red, green and 
blue, for generation of intensity patterns of arbitrary colours. 

30. A phase contrast imaging system (1) for synthesizing 
an intensity pattern I(x',y') of an image, comprising 

a source (2, 21, 41) of electromagnetic radiation for 

emission of electromagnetic radiation, 
a spatial phase mask (4, 23, 43) for phase modulation of 

electromagnetic radiation and having 

a plurality of individual resolution elements (x,y), each 
resolution element (x,y) modulating the phase of 
electromagnetic radiation incident upon it with a 
predetermined phasor value e*****^, and being 

positioned on a propagation axis of the electromagnetic 
radiation, 

means (5, 26, 44) for Fourier or Fresnel transforming the 
phase modulated electromagnetic radiation positioned 
on a propagation axis of the phase modulated radiation, 

a spatial phase filter (6, 27, 45) for phase shifting in a 
region of spatial frequencies comprising DC in the 
Fourier or Fresnel plane, the transformed electromag- 
netic radiation by a predetermined phase shift value 0 
in relation to the remaining part of the transformed 
electromagnetic radiation, 

means (7, 10, 26, 30, 44, 47) for forming the intensity 
pattern by Fourier or Fresnel transforming, 
respectively, the phase shifted Fourier or Fresnel trans- 
formed modulated electromagnetic radiation, whereby 
each resolution element (x,y) of the phase mask (4, 23, 
43) is imaged on a corresponding resolution element 
(x',y') of the image, 

the phasor values e 1 *^ of the phase mask (4, 23, 43) and 35 
the phase shift value 6 substantially fulfilling that 

for selected phase shift values 0, a being the average of the 4Q 
phasors e"*^ of the resolution elements of the phase mask 
(4, 23, 43), 

subject to the proviso that 

if 9«rt, the phase mask (4, 23, 43) is not divided into a 
matrix of rows and columns of resolution elements of 45 
the same size and shape, every fourth resolution ele- 
ment having the phasor value e* 5 and being distributed 
periodically and regularly across the area of the phase 
mask in such a way that every second row and every 
second column do not contain a resolution element with 50 
the phasor value e" 1 , the remaining resolution elements 
having the phasor value e'°, or, 

if 0«ji/2, the phase mask (4, 23, 43) is not divided into a 
matrix of rows or columns of the same size and shape, 
every second row or column having the phasor value 55 



means for selecting, for each resolution element, one of two 
phasor values which represent a particular grey level, and 

means for supplying the calculated phasor values e****^ to 
the elements (x,y) of the phase mask (4, 23, 43). 

32. A system (1) according to claim 30, wherein the 
intensity is zero for at least one resolution element (x 0 ',y^ 
of the intensity pattern, and wherein the phasor values e l *( x > y > 
of the phase mask (4, 23, 43) substantially fulfil that 



2kin- 



/(jc\ y ) = 2^1 ¥ mi^fa - y> + ?)] 



for selected phase shift values 0, fa being the phase of a. 

33. Asystem (1) according to claim 31, wherein the means 
for calculating the phasor values is adapted to calculate the 
phasor values e'*^ of the phase mask (4, 23, 43) in 
accordance with 



25 



30 



2sin- 



/(y,y) = 2[iTsiii(^-fly l y)+2]] 



for selected phase shift values 6, <fe being the phase of a. 

34. Asystem (1) according to claim 32, wherein the phase 
shift 0 is substantially equal to Jt and a is substantially equal 
to l A, and the phases <(>(x,y) substantially fulfil that 

/(y,y) = 2[i-cos(fly,y)j] 

J* Jsmtyi** y))dxdy = Q. 

Ax&y 



35. Asystem (1) according to claim 33, wherein the means 
for calculating the phasor values is adapted to calculate the 
phasor values e'*^ of the phase mask (4, 23, 43) in 
accordance with 

/(y, y) = 2[i- costly, yj)] 

JJsintylx, y))dxd y = 0. 

tixAy 

36. Asystem (1) according to claim 30, further compris- 
ing 

means for moving the region of spatial frequencies com- 
prising DC to a second part of the Fourier or Fresnel 
plane, and wherein 

the phase filter (6, 27, 45) is positioned at the second part 
of the Fourier or Fresnel plane for phase shifting of the 
transformed modulated electromagnetic radiation at the 
second part of the Fourier or Fresnel plane by 0 in 
relation to the remaining part of the electromagnetic 
radiation. 

37. Asystem (1) according to claim 36, wherein the means 
for moving the region of spatial frequencies comprising DC 
to a second part of the Fourier or Fresnel plane comprises an 
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optical component, such as a grating, a prism, etc, with an 
appropriate carrier frequency. 

38. A system (1) according to claim 36, wherein the means 
for moving the region of spatial frequencies comprising DC 

to a second part of the Fourier or Fresnel plane comprises the 5 
phase mask (4, 23, 43) in which the function of an optical 
component, such as a grating, a prism, etc, with an appro- 
priate carrier frequency, has been encoded. 

39. A system (1) according to claim 30, wherein the 
modulus of the Fourier transform of the phasors e 1 *^ at 10 
specific spatial frequencies have been adjusted to keep the 
intensity levels of the synthesized intensity pattern within a 
desired range. 

40. A system (1) according to claim 39, wherein the 
modulus of the Fourier transform of the phasors e'* c * ,y) at is 
specific spatial frequencies have been adjusted according to 

at least one of the following measures: 

a) adjusting the individual phasors e 1 *^-* of the resolution 
elements of the phase mask (4, 23, 43) maintaining 
prescribed relative intensity levels between intensities 20 
of resolution elements of the intensity pattern, 

b) adjusting the individual phasors e'*^ of the resolution 
elements of the phase mask (4, 23, 43) by histogram 
techniques, 

c) spatially scaling the phasor e***^ pattern of the phase 
mask (4, 23, 43), and 

d) utilizing half tone coding techniques. 

41. A system (1) according to claim 30, further compris- 
ing means for controlling power of the electromagnetic 30 
radiation in response to the intensity range of the intensity 
pattern. 

42. A system (1) according to claim 30, wherein each 
phasor e l<K * ,y) of the phase mask (4, 23, 43) is substantially 
equal to a selected phasor that has been selected from a set 35 
of two phasors with complementary phase values e^ 1 ^ 
and e 1 * 2 ^ in such a way that a specific spatial frequency 
distribution of the intensity of the electromagnetic radiation 

in the Fourier or Fresnel plane is attained. 

43. A system (1) according to claim 42, wherein the phase 40 
(ji(x,y) of phasors e 1 " 4 ^ of adjacent resolution elements 
alternates between the two possible complementary phasor 
values e* 1( ^> and c^ 2(x ' y \ 

44. A system (1) according to claim 30, wherein the phase 
filter (6, 27, 45) is shaped to match the spatial frequency 45 
content of the phasors e**^ of the spatial phase mask (4, 
23, 43). 

45. A system (1) according to claim 30, wherein modulus 
of the average value |a| of the phasors e'* ( ^ ranges from 0.1 

to 0.9. 50 

46. A system (1) according to claim 45, wherein modulus 
of the average value |a| of the phasors e***^ ranges from 
0.25 to 0.75. 

47. A system (1) according to claim 45, wherein modulus 

of the average value |a| of the phasors e**^ ranges from 0.4 55 
to 0.6. 

48. A system (1) according to claim 45, wherein modulus 
of the average value |a| of the phasors e l4,( *' y) is approxi- 
mately 0.5. 

49. A system (1) according to claim 45, wherein the phase 
shift 8 ranges from n/4 to 7ji/4. 



50. A system (1) according to claim 45, wherein the phase 
shift 6 ranges from to 3n/2. 

51. A system (1) according to claim 45, wherein the phase 
shift 8 ranges from 3x/4 to 5ji/4. 

52. A system (1) according to claim 30, wherein the phase 
shift 6 is approximately x. 

53. A system (1) according to claim 30, further compris- 
ing zooming means (10, 30, 47) for scaling of the intensity 
pattern. 

54. A system (1) according Lo claim 30, wherein the 
phase filter (6, 27, 45) comprises a spatial light modulator. 

55. A system (1) according to claim 30, wherein the phase 
mask (4, 23, 43) is adapted to perform the optical function 
of a Fourier-transforming lens by appropriate encoding of 
the phasors e 1 *^ of the phase mask (4, 23, 43). 

56. A system (1) according to claim 30, wherein the 
spatial phase filter (6, 27, 45) is adapted to perform the 
optical function of an output lens by appropriate encoding of 
the phase filter (6, 27, 45). 

57. A system (1) according to claim 30, wherein the 
source (2, 21, 41) of electromagnetic radiation is adapted to 
radiate electromagnetic radiation of different wavelengths 
corresponding to three different colours, such as red, green 
and blue, for generation of intensity patterns of arbitrary 
colours. 

58. A system (1) according to claim 30, further compris- 
ing a first and a second Fourier transforming lens (5, 7), the 
spatial phase mask (4, 23, 43) being positioned in the front 
focal plane of the first lens (5), the spatial phase filter (6, 27, 
45) being positioned at the back focal plane of the first lens 
(5), and the second lens (7) being positioned so that its front 
focal plane is positioned at the position of the back focal 
plane of the first lens (5). 

59. A system (1) according to claim 30, further compris- 
ing one Fourier transforming lens (44), the spatial phase 
filter (45) being positioned at the back focal plane of the lens 
(44). 

60. A system (1) according to claim 30, further compris- 
ing one imaging lens, the spatial phase filter (6, 27, 45) being 
positioned in the back focal plane of the lens. 

61. A system (1) according to claim 30, further compris- 
ing a polarizing beam splitter (24) and a quarter wave plate 
(25) and/or a phase filter (27) reflecting electromagnetic 
radiation incident upon it. 

62. A system (1) according to claim 30, wherein the 
spatial phase filter (6, 27, 45) changes the phase of the 
radiation in the region of spatial frequencies comprising DC 
and leaves the phase of the remaining part of the radiation 
unchanged. 

63. A system (1) according to claim 30, wherein the 
spatial phase filter (6, 27, 45) do not change the phase of the 
radiation in the region of spatial frequencies comprising DC 
and changes the phase of the remaining part of the radiation. 

64. A system (1) according to claim 30, wherein the 
spatial phase filter (6, 27, 45) blocks the radiation at the 
region of spatial frequencies comprising DC and leaves the 
remaining part of the radiation unchanged. 

65. A system (1) according to claim 30, wherein the 
source (2, 21, 41) of electromagnetic radiation is a Laser (2, 
21, 41). 
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It is certified that error appears in the above-identified patent and that said Letters Patent is hereby 
corrected as shown below: 



(Column 2. line 37), change "lexpIKjt^+jDl 2 " to -|exp[j(n/2)+j0]| 2 -. 
(Column 2. line 40), change "|expO(3n/2)+j0| 2 " to Hexp[j(3n/2)+j0]| 2 -. 
(Column 3, line 5) change "Zemike" to -Zernik-. 
(Column 6. line 59) change "e' a M -a" to -e w,y) -a - 
(Column 9. line 54) change "if to -1f-. 
(Column 10, line 59) change "if to -1f-. 

(Column 1 1 , line 60) change "where f j =Jf x 2 +// " to -where f=^f t 1 +// 

(Column 1 3. line 39) change " I o(x,y ) 1 12 =2[1-cos(0(x',y))]" to -I o(x'.y') 1 2 =2[1-cos(0(x',y))]-. 

(Column 13, line 41) change "phaseointensity" to -phase->intensity-. 

(Column 15, line 37) delete T'. 
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